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ABSTRACT 

We present a uniform X-ray spectral analysis of eight type-1 active galactic 
nuclei (AGN) that have been previously observed with relativistically broadened 
iron emission lines. Utilizing data from the XMM-Newton European Photon 
Imaging Camera (EPIC-pn) we carefully model the spectral continuum, taking 
complex intrinsic absorption and emission into account. We then proceed to 
model the broad Fe Ka feature in each source with two different accretion disk 
emission line codes, as well as a self-consistent, ionized accretion disk spectrum 
convolved with relativistic smearing from the inner disk. Comparing the results, 
we show that relativistic blurring of the disk emission is required to explain the 
spectrum in most sources, even when one models the full reflection spectrum 
from the photoionized disk. 

Subject headings: accretion, accretion disks - black hole physics - galaxies:nuclei 
- galaxies: Seyfert - X-rays:galaxies 

1. Introduction 

Black holes (BHs) are among the simplest objects in nature, able to be defined com- 
pletely by mass, spin and charge alone. But whereas the mass of a BH is relatively straight- 
forward to measure, provided material orbiting the BH can be observed (i.e., we are not 
blinded by the central engine), spin has proven much more challenging to constrain. With 
the advent of sensitive X-ray spectroscopy, allowing us to characterize processes very close 
to the BH in detail, we are now in a position to search for signatures of BH spin. This paper 
and its forthcoming counterpart (part II) represents the first survey to quantify the angular 
momenta of BHs in a sample of AGN using relativistically broadened iron lines as our spin 
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diagnostic. Here, in Part I, we examine the statistical need for including relativistic effects 
when modeling the emission from the inner accretion disks around BHs. 

As all current observables from a BH system are in one way or another related to the 
accretion process, any diagnostic of BH spin must be based on spin-induced changes to the 
accretion disk structure. BH spin manifests itself by setting the location of the innermost 
stable circular orbit (ISCO) which, in turn, sets the effective inner edge of a geometrically- 
thin accretion disk. For the purposes of this study, the "inner edge" of the disk is defined 
as the location inside of which no X-ray reflection signatures (e.g., iron line emission and 
Compton reflection hump) are produced either due to the disk becoming optically-thin or too 
highly ionized. As the magnitude of the (prograde) spin increases, frame-dragging stabilizes 
otherwise unstable orbits, causing the ISCO to move closer to the BH. The location of the 
inner edge of the disk can therefore be used to measure BH spin (e.g.. [Reynolds fc Nowak 
2OO3I : iBrenneman fc Revnoldslbood : iRevnolds fc Fabian! [20081 ). 



The position of the ISCO can be constrained through detailed modeling of the X-ray 
spectrum of the disk. Two popular modeling schemes include fltting the thermal continuum 
of the disk (e.g.. lNarayan et al. I l2007h and fitting the shape of the fluo rescent Fe K line emitted 



by th e disk in response to irradiation by the X-ray power-law (e.g.. IBrenneman fc Reynolds 
20061 ). The former method relies on a priori knowledge of the BH mass and distance as well 
as the disk inclination angle, however, and hence is most readily applied to Galactic BH 
systems (GBHs) in thermally dominant states as opposed to supermassive BHs found in the 
centers of AGN. By contrast, the shape of the iron line proflle depends upon the gravitational 
redshift and fluid velocity fleld, both of which depend only upon the mass and radius of the 
black hole (M/R). Because the radius of the event horizon is directly dependent of the mass 
of the BH, the iron line proflle is th erefore mass-independent, making it a useful technique 
for deriving spin in any BH system (iReynolds fc NowakI l2003l ) . This is an especially useful 



technique in AGN in which the BH mass is often poorly constrained. 

The shape of the iron line proflle is determined by several physical parameters: the BH 
angular momentum, the radial extent of the disk, the emissivity proflle of the disk and the 
inclination angle of the disk to the observer's line of sight. A proper model for the iron line 
must also take into account the effects of the Doppler shift. Special Relativity (e.g., beaming) 
and General Relativity (e.g., gravitational redshift), each self-consistently calculated for the 
given value of BH spin. The two widely available relativistic line models within the XSPEC 
spectral analysis software were computed only for single values of spi n; the diskline model 
describes only non-spinning Schwarzschild BHs (IFabian et al.lll989l ). and the laor model 
describes near-maximally (a = 0.998) Kerr BHs (lLaorlll99ll ). Thus, neither of these models 
allow for the self-consistent computation of a line proflle from an arbitrary spin BH, even if 
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one allows the inner edge of the line emitting region to vary. If one wishes to constrain BH 
spin using self-consistent models, therefore, it is necessary to create a new, fully-relativistic 
model code which enables the BH spin to be fit as a free parameter. 

We have developed such a model as both a line code (kerrdisk) and a convolution ker- 
nel (kerrconv) that can be applied to a full reflection spectrum expected from the diskQ To 
simplify the models, we have assumed that the BHs have no net charge (a very good approx- 
imation for all likely astrophysical BHs) and that there is a negligible amount of mass in the 
surrounding matter compared to the mass of the BH. Thus, the spac etime is described by the 
usual Kerr metric. These models were introduced and detailed in iBrenneman fc Reynolds 
(120061 ) . hereafter BR06. BR06 also contains detailed compariso ns of the kerrdisk and 



kerrconv models to similar codes developed by other groups (e.g.. iDovciak et al.l 120041 ). 



Our models assume that, in the co-moving frame of the disk surf ace, the ern itted 
line/reflection spectrum is limb-darkened using the same prescription as iLaon (119911 ). In 
fact, limb-darkening effects are usually completely unimportant to studies such as that con- 
ducted here; they only become significant for very high inclination angles, and furthermore 
it is only light-bending induced differences in the angle at which various parts of the disk 
are viewed that change the profile of the iron line. 

BR06 also presents a series of model fits to the 330 ks XMM-Newton/EPlC-pn spectrum 
of the canonical broad iron line Sy 1 AGN MCG-6-30-15, noting that the best fit is achieved 
with a kerrconv model acting on an ionized disk reflection spectrum, where the inner edge 
of the disk was r^[^ < 1.62 rg. These data have allowed us to place formal constraints on the 
BH spin in MCG-6-30-15 of a > 0.987 (90% confidence level). Because this fit was performed 
under the assumption of a strict truncation of the iron line region at the innermost stable 
circular orbit (ISCO) in the disk, this constraint is somewh at weakened (a > 0.92) if the iron 
line emission region does in fact bleed inside of the ISCO (IReynolds fc Fabian! |2008| ). Given 
that MCG-6-30-15 harbors the broadest iron line in an AGN observed to date, this model fit 
reinforces our expectation that this AGN contains a BH of near-maximal s pin. Many other, 
similar Sy 1 galaxies are also observ ed to have broad iron lines, however (iGuainazzi et al. 
20061 : iNandra et al.ll2007l : lMillerll2007l ). which begs two important questions: (1) can spectral 
observations of broad iron lines be used to robustly constrain BH spin, and (2) how is spin 
distributed between and among different types of BH systems? 

To begin addressing these questions, we have collected a sample of eight other Sy 1 AGN 
from the XMM-Newton archive that have previously been observed to have broad iron lines. 



^ These models are now available for public use in XSPEC12. Future updates will be made available initially 
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Using the same technique we used to model the spectrum of MCG-6-30-15 by first addressing 
the continuum and absorption and then considering the broad hne region, we have assessed 
each AGN independently, with the initial goal of establishing whether we can definitively 
show that relativistic effects from the inner accretion disk are evident in these spectra (Part 
I, this paper). Following this conclusion, we will examine the archived Suzaku spectra of 
the AGN in our sample in order to take advantage of the increased spectral resolution and 
extended energy range that this observatory offers. The goal of this second study is to obtain 
statistically robust BH spin constraints whenever possible (Part II, forthcoming). 

Section 2 of this work discusses the observations of the sources in our sample and the 
method used to reduce the EPIC-pn data and extract spectra. A summary of our modeling 
technique is given in §3. Section 4 contains the detailed spectral fitting results for each 
source. Results and tables comparing the relevant best-fit parameters for each source can 
be found in §5. Implications of these results, conclusions we have drawn and a prognosis 
for future work in this field are discussed in §6. Part II will discuss the Suzaku results and 
provide the first formal constraints on BH spin in our sample of AGN. 



2. Observations and Data Reduction 



Our source selection was based on work by lNandra et al.l (120071 ) and the review by lMiller 



( 120071 ). both of whom consider the robustness of broad iron line detections across different 
source populations. The Miller review, in particular, collects all the recent X-ray data 
published on Seyfert sources and categorizes these objects by the strength and robustness 
of their broad iron line features. We have used these two studies as a guide for selecting the 
candidates for our sample, which we list in Table [TJ 

For this work, our data have been collected from the XMM-Newton Science Archive, 
using only publicly available observations. XMM-Newton data were preferentially used as 
opposed to Chandra data because of the superior collecting area of XMM-Newton in the 
2 — 10 keV band where the broad Fe Ka feature is the most prominent. There were insufficient 
numbers of publicly available Suzaku datasets at the time Part I of this study began (early 
2007) to allow this work to be undertaken using Suzaku. The Suzaku study will be the focus 
of Part II in a forthcoming paper. 

In this paper, we exclusively focus on the EPIC-pn data to avoid cross-calibration issues 
between the pn, MOS and RGS instruments. Though the RGS spectra in particular can 
provide detailed insight into the 0.3 — 2.0 keV spectrum and help illuminate the nature of 
the soft absorption and excess seen in so many Sy 1 AGN, the mismatch in resolution and 
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signal-to-noise as well as the possibility of cross-calibration errors severely reduce the utility 
of including RGS data in our formal spectral fitting. The basic parameters characterizing 
the observations of each of our sources are shown in Table [H 

For each source in question we performed reprocessing (when necessary) and data reduc- 
tion with the SAS version 7.0.0 software, including the latest CCF calibration files. We have 
used the SAS epatplot task to compute the fraction of single, double, triple and quadruple 
events as a function of energy and compared these fractions to their nominal values as mea- 
sured from weak source observations. For sources that are affected by pile-up, these fractions 
deviate from the nominal values due to the higher probability of wrong pattern classification. 
No significant deviation from the nominal single and double distributions was found for any 
of our sources, indicating that our EPIC-pn observations are not affected by pile-up. Source 
and background spectra and light curves were extracted using the xmmselect task from the 
SAS GUI. Response matrices and ancillary response files were generated using the rmf gen 
and arf gen tasks. These were then grouped with the spectral files using the ftools package 
grppha with a minimum of 25 counts/bin in order to use a meaningful goodness-of-fit 

statistic. Spectral analysis was carried out as before for MCG-6-30-15 (BR06) with XSPEC 
version 12.4.0. 



3. Spectral Analysis Methodology 

In the following Section we outline the method used for the data reduction and spectral 
analysis of the AGN in our sample. We then discuss the spectral modeling of each source, 
culminating in a relativistically smeared ionized disk reflection spectrum model from which 
we infer the importance of taking GR into account when fitting each source. This is a critical 
first step toward measuring BH spin. 

3.1. Continuum Fitting 

We have endeavored to make the spectral analysis method as uniform as possible for 
all the AGN examined. There are, however, certain intrinsic physical differences between 
the sources that must be taken into account in order to properly model the continuum and 
isolate the iron line(s) in each case. Though all of our sources are Sy 1 AGN (or otherwise 
observed in states where they manifest broad disk lines), each system is unique in its physical 
properties. Some have higher signal-to- noise, either because of higher flux, longer observation 
time or both. Some exhibit evidence for complex, multi-zone warm absorbers (WAs) intrinsic 



AGN 


Date Obs. 


Exposure (s) 


Counts 


pn Mode 


Filter 


RA (h) 


Dec. (°) 


z 


Source Type 


MCG-5-23-16 


12/8/2005 


121578 


1.2 X 10^ 


LW 


Medium 


9.7945 


-30.7489 


0.0085 


Sy 1.9 


Mrk 766 


5/20/2001 


111789 


2.4 X 10^ 


PSW 


Medium 


12.3074 


29.8128 


0.0129 


Sy 1.5, NLSl 


NGC 3783 


12/19/2001 


131025 


1.6 X 10^ 


PSW 


Medium 


11.6505 


-37.7385 


0.0097 


Sy 1 


NGC 4051 


5/16/2001 


121889 


2.5 X 10^ 


PSW 


Medium 


12.0527 


44.5313 


0.0023 


Sy 1.5, NLSl 


Fairall 9 


7/5/2000 


29209 


2.0 X 10^ 


PFW 


Medium 


1.3961 


-58.8058 


0.0470 


Syl 


Ark 120 


8/24/2003 


112130 


2.3 X 10^ 


PSW 


Thinl 


5.2699 


-0.1502 


0.0327 


Syl T 


NGC 2992 


5/19/2003 


28917 


3.4 X 10^ 


PFW 


Medium 


9.7617 


-14.3264 


0.0077 


Sy 1.9 


3C 273 


6/30/2004 


19861 


3.2 X 10^ 


PSW 


Thinl 


12.4852 


2.0524 


0.1583 


QSO, Sy 1 



Table 1: Observation parameters for the eight AGN in our sample of broad iron line sources. Exposure times are after 
filtering, and counts are measured from 2 — 10 keV after filtering. 
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to the central engine, whereas some show only cold absorption from neutral hydrogen along 
our line of sight. A soft excess is seen in some sources but not others. And of course, the 
strength and breadth of the Fe Ka line varies from source to source as well, though all have 
been chosen because they have previously shown evidence for broad iron lines. 

In order to examine the iron line region in detail the continuum must first be accurately 
modeled. After excluding the energy ranges relevant to the iron line (4.0 — 8.0 keV) and 
the mirror edges (~ 1.5 — 2.5 keV), we examine the rest of the 0.6 — 10.0 keV spectrum, 
fitting it with a power-law continuum typical of an AGN which is then modified by cold 
photoabsorption from neutral hydrogen. We set the minimum value of A^h equal to the 
Galactic column density along our line of sight to the sourc^. If further absorption is 
necessary, we employ a second neutral absorber (either zphabs or zpcf abs) and allow the 
column (and covering fraction, if applicable) to vary as necessary to accommodate the cold 
absorption intrinsic to the system. Following the prior published analyses of these AGN, 
some sources require a second power-law component as well to properly model the basic 
continuum shape. We disregard energies below ~ 0.6 keV due to calibration uncertainties in 
this range for the EPIC-pn instrument. 

If significant residual features remain after fitting this continuum model which indicate 
the presence of a soft excess and/or warm absorption, we include these components one by 
one as long as they make a significant difference in the global goo dness-of-fit accord ing to the 
F-test: Ax^ < —4 for each new parameter introduced into the fit (jBevington 1969 ). The soft 



excess is typically parameterized by a blackbody component representing thermal emission 
in the central engine, but may also be modeled by bremsstrahlung emission or Comptonized 
emission from a thermal disk if either of these more complex forms gives a significantly 
better reduction in x^- The warm abs orption is modeled usin g the same xsTAR-generated 



multiplicative table model described in lBrenneman et al.l (120071 ) and BR06 in the analyses of 



NGC 4593 and MCG-6-30-15, respectively. Solar abundances are assumed for all elements 
and the redshift of the warm absorber (s) is set to that of the AGN host galaxy. Some AGN 
show no evidence for a warm absorber, while others statistically require up to two physically 
separate WA components, each exhibiting a distinct column density (A^h) and ionization 
parameter We note that the EPIC-pn does not have sufficient spectral resolution to 
discern velocity differences between the warm absorber and the host galaxy; thus we do not 
need to include the warm absorber velocity as a free parameter. 
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See 



http://heasarc.Ksfc.nasa.KOv/cgi-bin/Tools/w3nh/w3nh.pL for the column density calculator as well as 
Dickey fc LockmanI (|l990f ) and lxalberla et al.l |2005h . 
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Once the continuum has been properly modeled, we initially freeze the parameter val- 
ues for all components except the power-law spectral index (F) and normalization (flux in 
phcm~^s~^) and restrict our energy range from 2.5 — 10.0 keV in order to focus on the 
hard X-ray spectrum and the iron line region. Including the energies from 4.0 — 8.0 keV 
again, we check for residual emission (or absorption) lines from 6.4 — 6.97 keV, indicating 
the presence of unmo deled neutral and/or partly ionized Fe Ka in the spectrum. Because 
these sources have been pre-determined to possess significant, broadened iron emission, it 
is not surprising that such residual emission features are seen in each case. We begin by 
attempting to fit the 6.4 keV line of neutral iron (and any other ionized iron lines) with a 
simple Gaussian feature, with the line width frozen at cr = OkeV (i.e., intrinsically narrow) 
and the rcdshift again set to the cosmological value for the host galaxy. Such narrow lines 
are expected to originate from the outer accretion disk, or perhaps from reprocessing in a 
Compton-thick toroidal region surrounding the nucleus. As expected, in each case we note 
significant residual "wings" remaining around the Fc Ka feature after this narrow core was 
fit, indicating the presence of a broad component to the line originating from closer to the 
BH. 



3.2. Fitting the Iron Line Profile 

To assess the robustness of the broad iron component and the importance of relativistic 
smearing in the spectrum, we follow our continuum fitting procedure (including modeling 
the narrow core of the iron line) by running each source through two additional steps of more 
complex modeling in the broad iron line region. 

We begin by adding to the spectral model a broad iron line described by: (1) a broad 
Gaussian line (a =free), replaced by (2) the diskline model (spin hardwired at a = 0.0), 
then replaced by (3) the laor model (spin hardwired at a = 0.998). The broad Gaussian is 
used primarily to gauge the energy and velocity width of the line in order to infer its radius 
of origin, whereas the relativistic line models are thought to be more physically accurate 
representations of the fine morphology. The energy of the line is initially set to a rest- 
frame energy of 6.4 keV, corresponding to the energy of the cold Fe Ka line. We use the 
change in x^/^ with the addition of the broad line component to verify whether a broad 
line is statistically warranted in our global model. According to this criterion, all of our 
sources required a broad line, and relativistic broadening was robustly preferred over a 
simple Gaussian in half of the cases. This is unsurprising since our objects were preselected 
to possess reported broad iron lines. Using the simple diskline and laor models, we also 
note the inner radius of disk emission: a broad line with the potential to diagnose BH spin 
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should have an inner disk radius of roughly rmin < 6 r^. We relax this restriction here to 
allow for differences in fit that may be achieved through different models, requiring that 
^min ^ 20 Tg to infer that relativistic broadening is important. 

We next perform a study to assess the degree to which the X-ray reflection continuum 
may be contaminating (or even mimicking) the broad iron line. To begin w 'ith, we replace 



the b road iron line with a static (unsmeared) ionized reflection spectrum (IRoss fc Fabian 



20051 ) and refit. We then note the change in the goodness-of-fit criterion when the reflec- 
tion spectrum is subjected to smearing characterizing a relativistic disk by convolving the 
spectrum with a line profile smearing kernel. The first such kernel we use is for a non- 
spinning BH (a = 0.0) and is known as rdblur. The second convolution model assumes a 
maximally-spinning BH (a = 0.998) and is known as kdblur. 

We force the outer radius of disk emission to be r^ax = 400 rg in all the relativistic 
broad line fits. The value itself is somewhat arbitrary; as long as a > 2 (where a here 
denotes the emissivity index of the disk and the disk radiates as a function of radius: r~°), 
emission from this outer part of the disk is negligible, and we have already incorporated a 
separate narrow line core component whose width dictates an origin in the outer disk or 
torus region. For the ionized disk component, we begin with the assumptions of solar iron 
abundance and relative neutrality (Fe/solar = 1 and C, = 30ergcm~^ s^^ ). We then relax 
these assumptions and allow these two parameters to fit freely, provided that reasonable 
statistical error bars can still be obtained for each parameter. We also note that all the 
models used above assume that the disk radiation is truncated at the ISCO; that is, that no 
contribution to the continuum or iron line profile originates from within this region. 



4. Source-by-source Model Fitting 
4.1. Target selection 



As mentioned in the previous Section, we have selected the AGN for our sample pri- 
marily based on the research performed by other groups on the robustness of broad iron line 
observations in various AGN and the likelihood of bei ng able to identify r e lativistic blu rring 
effects in t hese sources using the iron line method (INandra et al.l 120071 : iMilleii 120071 ). In 
particular, iMilleii (120071 ) discusses recent results from 30 Seyfert AGN in which relativistic 
disk lines have been detected and published, and divides those sources up into three "tiers" 
based on the nature and robustness of the detections. Many of the sources overlap with 
those presented in Nandra et al. Though we have examined observed spectra from all of the 
9 AGN listed in Tier 1 (the most robust cases), NGC 3516 and NGC 4151 were not included 
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in ou r study due to their extremely complex soft spectra (ITurner et al.ll2005l : ISchurch et al. 



2004h . We did analyze the other 7 sources in full: 3C 120, MCG-6-30-15, MCG-5-23-16, 



NGC 2992, NGC 4051, NGC 3783 and Mrk 766. 3C 120 is not included in this work, how- 
ever, as our spectral fitting results for this object were inconclusive in several ways due in 
part to a paucity of photons. We have also examined 3 other sources from Tier 3 that are 
mentioned prominently in Nandra et al. : Fairall 9, Ark 120 and 3C 273. 



4.2. Source-by-source results 

We present the best-fitting results of our spectral analysis for each source in its own 
subsection below. In many cases the kdblur and rdblur results were very similar, implying 
that little can be inferred about the spin of the BH in these AGN without more data. A 
more detailed comparison of all of our models is taken up in §51 Note that the reflection 
fractions quoted are not fitted values, but rather are estimates calculated for the purpose of 
comparison to other studies which may model disk reflection differently (e.g., with pexrav 
rather than ref lion). We expand on this point in greater detail in J|5l Detailed spectral 
fitting results for MCG-6-30-15 can be found in BR06 and are therefore not included herein. 

A table of best-fit parameters for each source is included, as are figures representing the 
residual features left behind after initially fitting the spectrum with an absorbed power-law, 
the final best fit to the spectrum, and the best-fit model components for each object. 



4.2.1. MCG-5-23-16 



MCG-5-23-16 is a moderately absorbed Seyfert galaxy of intermediate type (Sy 1.9) 
It is relatively nearby at a redshift of z = 0.0085, and with a typical 2 — l OkeV flux of 
B X 10^^^ erg cm~^ s~^ it is one of the brightest known Seyfert galaxies (iReeves et al. 



2-in 



20061 ) . The source has been observed pre viously to have a soft excess below ~ 1 keV and 



an absorbing column of Nn ~ 10^^ cm ^ (jPewangan et al.ll2003l : iBalestra et al.ll2004j ). and 



ASCA observations have indicated the presence of a broad Fe Ka line of EW ~ 200 eV 
(IWeaver et al.lll997l . ll998l ). The line was successfully modeled with a narrow core at the rest- 
frame energy of 6.4 keV and a broad component superposed on it. This feature was modeled 
with similar success in the XMM-Newton observations of Dewanger i et al. and Reeves et al. , 
as well as the Suzaku observation also published by these authors (IReeves et al.l 120071 ). 



In December 2005, MCG-5-23-16 was simultaneously observed with XMM-Newton^ 
Chandra^ Suzaku and RXTE. The XMM-Newton results are reported by Reeves et al. as well 
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as Br aito et al. , and our data reduction followed those works (iReeves et al.ll2006l : iBraito et al. 



20061 ). The EPIC-pn instrument had a net exposure of 96 ks and returned ~ 1.2 x 10 photons 



after filtering. 

As suggested by Reeves et al. , a simple photoabsorbed power-law fit does not model the 
continuum well, especially below ~ 1 keV. Adding in a source of soft thermal emission also 
leaves prominent residuals and does not adequately account for the shape of t he soft excess. 



so fol lowing the lead of Reeves et al. we have employed a two-power-law model (IReeves et al. 



20061 ): one of the power-law components is affected only by Galactic photoabsorption (iVn = 
8 X 10^° cm~^ as per Reeves et al. ), and hence leaves the AGN system effectively unscattered, 
while the other component is also subject to absorption intrinsic to the AGN, in this case 
with a fitted absorbing column density of iVn = 1.19±0.01 x 10^^ cm~^ . This component also 
experiences scattering within the system. We set both power-law photon indices to be equal, 
indicating that the two power-law components originate from the same physical reservoir of 
photons (with F = 1.66 ± 0.01). Calculating the ratio of the flux of the scattered power-law 
component to the unscattered power-law component (from 0.6 — 10.0 keV, fitting only the 
continuum) yielded an estimate of the optical depth of the scattering plasma. In this case 
T = 7.37 X 10"'^ ph cm~^ s~^ /1. 98 x 10~^ ph cm~^ s~^ ~ 0.04, so the scattered fraction is low, 
implying that the electron plasma is optically thin. No evidence for spectral curvature due 
to the presence of a WA was found in this source. 

The spectrum does show evidence for two prominent, narrow emission lines of iron Ka 
and Kp (at = 6.40 ± 0.01 keV and EKf3 = 7.06^°;°? keV, with EWxa = 167l}^eV and 
EWxfs = 104^30 eV, respectively), as well as three narrow absorption lines of intermediately 
ionized iron or nickel at Ei = 7.241°;°^ keV, E2 = 7.48 ± 0.04 keV and E3 = T.Sbtof^keV 
{EWi = -79+^^ eV, EW2 = -UlttleV, and EW^ = -143 ± 43 eV, respectively). 

After successfully modeling the continuum and narrow line parameters for MCG-5-23- 
16, we restricted our attention to the hard spectrum (2.5 — 10.0 keV) and analyzed the broad 
line component of this source using the method outlined above in §3.21 The initial broad 
Gaussian fit implied that the broad line originates in neutral iron with E = 6.34 ± 0.02 keV 
and a = 0.28 ± 0.03 keV. The width of the line suggests an origin within d ~ 100 rg of 
the BH. A laor line component fit the data marginally better than a diskline, with an 
improvement of Ax^/Az/ = —145/ — 4 in global goodness-of-fit over a model without a broad 
iron line component. Here, u is the number of degrees of freedom in the fit, defined by the 
total number of spectral channels in the data minus the number of free model parameters, 
so a decrease in u implies an increase in the number of free model parameters. 

We determined that an ionized disk reflection spectrum convolved with kdblur provides 
the best model fit, though the inner disk radius could only be weakly constrained (rmin < 
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63 Tg). This prevents us from being able to place any meaningful constraints on the spin of 
the BH in this source with these data, in spite of the presence of relativistic broadening in 
the iron line profile. Note that eliminating soft energies from the fit removed the statistical 
need for a second, soft power-law in the fit. For a full listing of best-fit model parameters and 
error bars for MCG-5-23-16, see Table [2l The residual iron line feature, best-fit model and 
best-fit model components for the hard spectrum are shown in Fig. [H The diskline and 
laor fits are shown in Table [TOl while comparisons of the ionized disk reflection spectrum 
with and without relativistic smearing can be found in Table [TTl 



4.2.2. NGC 3783 



NGC 3783 is a bright, nearby Sy 1 galaxy at a redshift of z = .0097. It was first 
detected in X-rays with the Ariel- V all-sky survey (IMcHardv et al.lll981l) and subse quently 
in the high Galactic latitude survey conducted by HEAO-1 (jPiccinotti et al.lll982l ). Since 
these early detections, there have been many observations of NGC 3783 with higher res- 
olution instruments: A ROSAT observation of the source showed evidence of an ionized 



absorber in the soft band (ITurner et al.lll993l ). which was confirmed during ASCA obser- 
vations (iGeorge et al.lll995l . Il998l ). High-resolution grating observations of NGC 3783 with 
Chandra and XMM-Newton followed, unveiling the soft spectrum of this source in detail 
jKaspi et al.lboool . boOll . I2OO2I : Isiustin et aPbooi ISehar et al. l l2003h . The higher signal-to- 
noise of XMM- Newton, in particu lar, has also enabled the iron line to be studied extensively 
in this source (IReeves et al.ll2004l ). Using two observations from December 2001, Reeves et 
al.have obtained ~ 240 ks of data on NGC 3783. Their global fits to the merged EPIC-pn 
spe ctrum are currently second in length only to the ~ 330 ks observation of MCG-6-30- 
15 (IFabian et al.ll2002l ). In these observations, the source has an average 2 — lOkeV flux 
of = 6.8 X 10~^^ ergcm"^ s~^ . The spectrum is noted to have iron line peaks at 

6.39 ± 0.01 keV and 7.00 ± 0.02 keV, representing neutral Fe Ka {EW = 123 ± 6eV) and 
a blend of ionized Fe Ka and Fe Kj3 [EW = 34 ± 5eV), respectively. A strong absorption 
line of intermediately ionized iron is also observed at 6.67 ± 0.04 keV {EW = 17 ± 5eV) 
that exhibits variability in direct correlation with that of the continuum flux. A weak red 
wing of the 6.39 keV line is noted by the authors, but once warm absorption is taken into ac- 
count in the systern . the requirements for a broad Fe Ka component are significantly reduced 
jReeves et al.ll2004h . 



We analyzed the first of the two December 2001 XMM-Newton observations. We chose 
not to merge the event files of both data sets to avoid the uncertainties inherent in so doing, 
instead focusing on only one data set with an effective exposure time of ~ 120 ks and a 
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Fig. 1. — Top left: The 2.5 — 10.0 keV spectrum of MCG-5-23-16 fit with a two-power-law 
model modified by Galactic photoabsorption; data-to-model ratio. Note the prominent Fe 
K line and the large residual emission component above 7keV, which likely indicates the 
presence of ionized disk refiection. 'x^/v = 2675/1499(1.78). Top right: Best-fit model for 
MCG-5-23-16, including our continuum model, two narrow emission lines, and three narrow 
absorption lines as well as an ionized disk refiection spectrum convolved with our kdblur 
relativistic smearing kernel. xV'^ = 1485/1488(1.00). Bottom: An E^F{E) plot depicting 
the relative fiux in each of the best-fit model components in the spectrum. 
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Table 2: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of MCG-5-23-16, 
including components and parameter values from 0.6—1.5 keV for completeness. The energies 
from 1.5 — 2.5 keV were not included due to the presence of absorption edges from the XMM- 
Newton mirrors. Error bars are quoted at 90% confidence. For the zgauss lines, we required 
each to be intrinsically narrow, i.e., a = OkeV. Redshifts were frozen at the cosmological 
value for the source, in this case z = 0.0085. Note the absence of warm absorption and a soft 
excess in this object. When no error bars are quoted, the parameter in question is frozen at 
the given value. 
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total of ~ 1.6 X 10^ photons after filtering. This number of counts is still large enough to 
obtain valid statistical fits to the model parameters used. We follow the data reduction 
steps of Reeves et al. , excluding the last part of the observation due to contamination from 
a background flare. Our analysis follows in the same spirit as these authors, but in order to 
remain consistent with the method described above in ^we do not make an y assumptions 



about the soft spectrum and continuum emission based on prior RGS results (IBlustin et al. 



2002|). 



We began our analysis of the time-averaged spectrum with a simple photoabsorbed 
power-law fit using the Galactic absorbing column of A^h = 8.5 x 10^° cm~^ , and neglecting 
the energy ranges associated with the mirror edges and the iron lines, as in ||3l This fit 
leaves obvious residuals in all parts of the spectrum, however, and is especially poor for the 
soft energies below ~ 1.5 keV. Continuum curvature associated with a soft excess and warm 
absorption are both clearly evident. As discussed by Reeves et al. , the fit was much improved 
with the addition of a thermal blackbody component {kT ~ 0.07±0.01 keV, as compared with 
their value: kT ~ 0.09 ± 0.01 keV). We approached the question of the warm absorption in 
the same manner as Reeves et al. , using our xstar table model to parameterize the column 
density and ionization level of the absorbing medium. As with these authors, we found 
statistical evidence for a multi-zone warm absorption structure, though our fitted values for 
the column densities and ionization parameters of these zones vary from those of Reeves et al. 
These authors find A^hi = 1-1 x 10^^ cm~^ and log^i = —0.1 (for this zone, Fe/solar = 10, 
in contrast to the other zones with solar iron abundance), A^h2 = 1.2 x 10^^ cm~^ and 
log ^2 = 2.1, with A^H3 = 4.4 x 10^^ cm~^ and log^s = 3.0 (no final error bars were given). 
By contrast, we find A^hi = l-SOl^^Jg x 10^° cm'^ and log^ = 0, A^h2 = l-^Stosl x 10^3 cm'^ 
and log 6 = 2.82 ±0.94, with A^hs = SA^tof^ x 10^3 cm'^ and log ^3 = 3.79+^4- As with the 
results of Reeves et al. the third, neutral absorber of smaller column density was required to 
obtain a quality fit, but a very high iron abundance was necessary: Fe/solar = 10, where the 
abundance was approximated iteratively and held constant in the final fit because it was at 
the endpoint of its parameter space. We used the same technique to constrain the ionization 
parameter of this absorbing zone for the same reason. The differences between our results 
and those of Reeves et al. are n ot surprising, howev er, because Reeves et al. base their fits on 



the RGS results for this source (IBlustin et al.ll2002l ). We do not use any a priori information 



to augment or guide the EPIC-pn spectral fits due to calibration uncertainties. Finally, we 
allowed for a second neutral hydrogen absorbing column intrinsic to the source (using the 
zpcf abs model component) in an effort to improve the goodness-of-fit and found that the 
spectrum preferred a higher value than the Galactic column: A^h = 2.77I0I3 ^ ^0^^ cm~^ with 
a covering fraction of / = 87±2%. This was also in contrast to the Reeves et al. result, where 
the cold hydrogen column was held at the Galactic value and no intrinsic neutral hydrogen 
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absorption was allowed. The differences in absorbing structures between our model and that 
of Reeves et al. also explain the difference in photon indices of the power-law continuum in 
the two models: our best-fit value was marginally softer, at F = 1.831q;o4, whereas Reeves 
et al. found a harder photon index of F = 1.73 ± 0.01. 

Reeves et al. found two emission peaks in the spectrum at 6.39 ± 0.01 keV {EW = 
123 ± 4eV) and 7.00 ± 0.02 keV {EW = 34 ± 5eV), representing cold Fe Ka and likely 
a blend of ionized Fe Ka and Fe K/3, respectively. The latter line, in particular, showed 
no appreciable variance over the observation, leading the authors to postulate a relatively 
distant origin for the lines away from the central parts of the accretion disk. A 6.67±0.04 keV 
absorption feature {EW = 17 ± 5eV) of highly ionized iron was also seen which did appear 
to vary with time and was strongest when the continuum flux was highest, suggesting an 
origin in the region of the warm absorber within 0.1 parsec of the nucleus. Though inclusion 
of the WA structure (tailored to include this absorption feature) did lessen the statistical 
case for a broad iron line in NGC 3783, Reeves et al. nonetheless identified a residual broad 
feature which they fit using a diskline model with a = 3.3 ± 0.5, i = 19 ± 9°, rmin = 6.0 rg 
(held constant, with rmax = 100 rg), and EW = 58 ± 12 eV. 

We also detected the emission and absorption features discussed by Reeves et al. , though 
we found that the equivalent widths for the Fe Ka and Fe Ka/KjS blend differed from the 
fits performed by these authors: EWxa = 65 ± 12 eV, EWko/k/s = 28 ± 3 eV. As mentioned 
above, the inclusion of the three- zone WA structure eliminated the need to include a separate 
Gaussian absorption line at 6.67 keV. It should also be noted that, contrary to the Reeves 
et al. analysis, we required the core of the 6.4 keV line to be intrinsically narrow, so it 
is not unusual that we obtain a smaller equivalent width than Reeves et al. When we 
included a broad diskline component of the neutral Ka line, the goodness-of-fit improved 
dramatically: A^^/Az/ = —170/ — 4. This fit proved superior to both a laor line model 
and a broad Gaussian component, the latter of which particularly did not model the shape 
of the line well. Our best fit to the hard spectrum of NGC 3783 convolved an ionized 
disk reflection spectrum with rdblur, as can be seen in Table [3] below. The inclusion of 
the ionized disk reflection spectrum component (reflion) accounted for the soft excess 
self-consistently, without the need to include the phenomenological blackbody component. 
Adding in the effects of relativistic smearing on this reflection spectrum, our constraints 
on the radial extent of the disk showed that r^in < 46 rg, making it difficult to distinguish 
a rapidly-spinning BH from a slowly-spinning one, even though relativistic disk effects are 
clearly important in broadening the iron line profile. The residual iron line feature, best- 
fit model and best-fit model components for the hard spectrum are shown in Fig. [21 The 
diskline and laor fits are shown in Table [TOl while comparisons of the ionized disk reflection 
spectrum with and without relativistic smearing can be found in Table [TTl 
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Fig. 2. — Top left: The 2.5 — 10.0 keV spectrum of NGC 3783 fit with a power-law model 
modified by Galactic photoabsorption, along with a three-zone warm absorber and a black- 
body soft excess. xV^^ = 2751/1487(1.85). Top right: The 2.5 - lO.OkeV best-fit model 
for NGC 3783, including our continuum model and two narrow emission lines as well as 
an ionized disk reflection spectrum convolved with the rdblur relativistic smearing kernel. 
X^/v = 1480/1487 (1.00). Bottom: An F(E) plot depicting the relative flux in each of the 
best-fit model components in the spectrum. 



-18- 



Model Component 


Parameter 


Value 


phabs 


Nn ( ) 


8.50 X 10^0 


zpcf abs 


AT { / >-v ~2 \ 

iVH2 [ cm j 

Jcov I /Oj 




2.77+°-l^ X 1021 
88^2 


WA 1 


^VwAi (, cm ) 

log CwAl 




1.80+°-23 X 10^0 



WA 2 


V WA2 (, cm j 

iOg ^WA2 




1.58+°-^^ X 10^3 

9 09+0.94 


WA 3 


AT ( n-rr. —2 \ 

^VwA3 [ cm j 

log^WAS 




8.43+°-°^ X 10^3 

70+0.02 

— U. 14 


zpo 


r 

po 




1 oQ+0.03 

X . 00 A r\A 
— U.U4 




flux ( ph cm-2 s" 




1.33^°:°^ X 10-4 


zgauss 


hj ke V j 








a / I, — 2 - 

nux ( pn cm s 


) 


3.2512-37 X 10-^ 




EW ( eV) 




651^ 


zgauss 


E ( keV) 




7.01+°°2 

' " — 0.03 




flux ( ph cm^^ 


-1 ^ 


l.ll+°-2o X 10-^ 

—0.20 




EW ( eV j 




28^5 


rdblur 


a 

^min (^g) 
^max ('"g) 




1.90+1-^° 

— U.o7 

23+0 

1 7+29 
i 1 

400 


ref lion 


Fe / solar 




0.83l°i2 




^refl (erg cm^^s" 




30+24 




Trefl 




1 Oq+0.03 

1.00_Q Q4 




Rred 




o.49;°:22 




flux ( ph cm~2 S" 




7.4612;^? X 10-8 






1480/1487(1.00) 



Table 3: Best-fltting model parameters for the 2.5—10.0 keV spectrum of NGC 3783, including 
components and parameter values from 0.6 — 1.5 keV for completeness. The energies from 
1.5 — 2.5 keV were not included due to the presence of absorption edges from the XMM- 
Newton mirrors. Error bars are quoted at 90% confidence. For the zgauss lines, we required 
each to be intrinsically narrow, i.e., a = OkeV. Redshifts were frozen at the cosmological 
value for the source, in this case z = 0.0097. Note the multi-zone warm absorber present in 
this object, similar to many other Sy 1 sources. 
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4.2.3. Mrk 766 



Mrk 766 is designated as a classic, bright NLSl galaxy (Sy 1.5) with a redshift of 



0.0129 and a typical 2 - lOkeV flux of F2-10 ~ 2.5 x lO'^^ ergcm-^ s~i flPounds et al. 



20031 ) . Previous X-ray observations of this source have provided contradictory evidence 
on the detection of a broad Fe Ka feature. Although not one of the most convincing cases, 
Mrk 766 was included in an ASCA spectral su rvey of bright Seyfert galaxies showing evidence 



of relativistic iron lines (INandra et al 



19971). A separate analysis of simultaneous ROSAT 



and ASCA observations (ILeighly et al.lll996l ). however, showed the X-ray spectrum to be 
described by a power-law of index increasing strongly with flux from F ~ 1.6, but with 
only a narrow Fe Ka emission line {EW ~ 100 eV at 6.4 keV). A later observation with 
BeppoSA X found a steeper power-law (F ~ 2.2), and evidence for an absorption edge at 



~ 7.4 keV (IMatt et al.ll2000l ). implying strong reflection from intermediately ionized material. 
Interestingly, based on XMM-Newton observations of the source. Turner et al. found that 
energy-time maps of Mrk 766 reveal a periodic energy shift in an ionized component of Fe 
Ka emission, with a period of ~ 165 ks. This can be interpreted as evidence for emission 
from orbiting gas within ~ 100 rg of the central BH. A likely explana tion is that this ga s 
represents a hot spot on the disk illuminated by magnetic reconnection (ITurner et al.ll2006l ). 



We examine the May 2001 XMM-Newton observati on of Mrk 766 take n by Mason et 
al. , from which the RGS results were published in 2003 jMason et al]l2003l) and the EPIC 



o.69i°:°D 



to moderately ionized {Nm = 2.73lJ];^2 x 10^^ cm 



spect ra were explored more thoroughly by other authors (iPounds et al.l 120031 : [Turner et al. 
20061 ). In our re-analysis, the EPIC-pn data have an effective exposure time of 128 ks and 
the filtered data yield ~ 2.4 x 10^ photons. 

Mrk 766 bears many spectral similarities to MCG-6-30-15, showing significant complex- 
ity beyond a simple photoabsorbed power-law fit. There is clear evidence for a thermal soft 
excess below ~ 1 keV, which we initially parameterize with a phenomenological blackbody 
of kT = 0.18 

-aoi keV. There is also evidence for two distinct physical zones of low (A^hi — 
3.29lJ]i^ X 10^1 cm-2, log^i 
2 = 2.17+o;2o) 

tion by the Milky Way (Ah = 1-71 x 10^° cm~^ ). We find that the power-law in this source 
has a spectral index of F = 2.13 ± 0.01, which is s lightl y harder than the value determined 
for the time-averaged spectrum by [Pounds et al.l (120031 ) (F = 2.21 ± 0.01). Unfortunately, 
these authors did not probe the spectrum below 3 keV, so we cannot compare our soft excess 
or warm absorption parameters to theirs. 

The narrow Fe Ka core of Mrk 766 has a rest-frame energy of 6.40 ±0.01 keV, consistent 
with neutral iron, and an EW = 47 ± 8 eV, consistent with the strength of the narrow line 
of neutral iron found by Pounds et al. Fitting this line with a Gaussian component leaves 



log ^2 = '^■^7-0^1) intrinsic absorption. These components exist in addition to cold absorp- 
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significant residuals strongly indicative of the presence of a broad line. When this residual 
feature is fit with a broad Gaussian component its rest-frame energy is indicative of ionized 
iron: E = 6.62 ± 0.07 keV. Its width (a = 0.30 ± 0.08 keV) is indicative of an origin 
~ 100 Tg from the black hole. This broad component is most successfully fit here with a 
kdblur (ref lion) model (removing the now- unnecessary blackbody component), yielding 
comparable results to the disk reflection model employed to model the broad iron line in 
Pounds et al. (reduced = 1.09 in our work, compared with 1.11 in Pounds et al. ), though 
with slightly tighter parameter constraints (e.g., rmin < 14 rg in our work vs. 70 ± 12 rg in 
Pounds et al. ). Our model also detects no absorption features above 8keV, as discussed 
by Pounds et al. , possibly due to our different approach toward modeling the reflection 
spectrum, i.e., the use of ref lion vs. the xion model employed by Pounds et al. We 
note, however, that there is no significant difference between our kdblur and rdblur fits, or 
between out diskline and laor fits, indicating that we are not able to distinguish between 
a spinning and non-spinning black hole in this source using these data. The final best-fit 
results are shown in Table HI Fig. [3] shows the iron line residual, best fit to this residual, and 
relative contributions of the individual model components, respectively. The diskline and 
laor fits are shown in Table (TUl while comparisons of the ionized disk reflection spectrum 
with and without relativistic smearing can be found in Table [Til 



4.2.4. 3C 273 



3C 273 is the most distant source in our sample at a redshift of z = 0.1583. Classified as 
a bright variable quasar, this AGN displays a strong jet and during epochs of radi o-loudness it 



posse sses the flat X-ray spectrum of a blazar, with highly beamed jet emission ( iTiirler et al. 



20061 ). When the jet is reduced in strength, however, this source has been observed to exhibit 
Sy 1-type accretion disk signatures such as a broad iron line. A jet-minimum state occurred 
in March 1986 and allowed Robson et al. (1986) to identify a new near-infrared spectral 
component. An even better opportunity arose in early 2004, when the sub-millimeter flux of 
3C 273 was observed to be almost two times lower than in 1986. This new minimum triggered 
a slew of simultaneous observations with instruments in all wavebands such as INTEGRAL, 
XMM-Newton and RXTE, among several other optical, radio and sub-millimeter telescopes. 
The 2 — 10 keV flux during this period was E2-10 = 6.7 x 10^^^ erg cm"^ s^^ . S uch a low flux 
in 3C 273 has only been measured twice in t he past, by Ginq a in July 1987 (ITurner et al. 



19901) and by BeppoSAX on 18 July 1996 flHaardt et al.l 119981 1 . the latter of which was 
coincident with the low sub-mm flux mentioned above. This X-ray/sub-mm correlation 
strongly supports a synchrotron self-Compton origin for the X-ray jet emission in 3C 273. 
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Fig. 3. — Toj) left: The 2.5 — 10.0 keV spectrum of Mrk 766 fit with a continuum composed of 
a power-law model modified by Galactic photoabsorption, a two- zone warm absorber model 
and a blackbody soft excess, y^^ jv = 1681/1315 (1.28). Top right: The 2.5 — 10.0 keV best-fit 
model for Mrk 766, including our continuum model and an ionized disk reflection spectrum 
convolved with the kdblur relativistic smearing kernel. xV'^ = 1426/1308 (1.09). Bottom: 
An F(E) plot depicting the relative flux in each of the model components in the spectrum. 
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Table 4: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of Mrk 766, includ- 
ing components and parameter values from 0.6 — 1.5 keV for completeness. The energies 
from 1.5 — 2.5 keV were not included due to the presence of absorption edges from the 
XMM-Newton mirrors. Error bars are quoted at 90% confidence. We required the zgauss 
line representing the core of the Fe Ka feature to be intrinsically narrow, i.e., a = OkeV. 
Redshifts were frozen at the cosmological value for the source, in this case z = 0.0130. A 
two-zone warm absorber is present in this object, as in MCG-6-30-15 and many other Sy 1 
sources. We were not able to constrain the iron abundance for the ionized disk model, so we 
elected to fix it at the solar value. 
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Tiirler et al. obtained a 20 ks XMM-Newton observation of 3C 273 during this jet- 
minimum state in June 2004. We use their (reprocessed) thin-filter EPlC-pn observation 
in the interest of collecting as many photons as possible, though due to photon pile-up it 
i s necessary to exc lude the centralmost region of the source, as detailed by the authors 
( iTiirler et al.ll2006l ). In our filtered data set we capture ~ 3.2 x 10^ photons. 



Tiirler et al. note the inadequacy of a simple photoabsorbed power-law fit to the data. 
Their best fit to the continuum is achieved using two power-law components: Fhard = 1-63 ± 
0.02 and Fgoft = 2.69 ± 0.06, with the hard component flux ~ 2.3 times the soft flux. Both 



components were modified by Galactic photoabsorption with A^h = 1-79 x 10 cm . We 
based our initial continuum fit on theirs, with two power-law components: Fhard = 1-72" 



-0.11 
-0.28 



and Fgoft = 3.01_q'7]^, both consistent with the Tiirler et al. results. The ratio we calculate 
between the fluxes of the hard and soft power-law components is ~ 4.18, however. This is 
significantly higher than that calculated by Tiirler et al. , and the error bars on the spectral 
indices, in particular, led us to consider ehminating the soft component. We did so, but 
still found it necessary to parameterize the remaining soft excess with a thermal bbody 
model at first (kT = 0.14to;o2 keV). Such a component is not entirely phenomenological in 
nature, however: in many radio-loud sources a soft thermal excess has been noted, and is 
often interpreted as ar ising from the interaction of the jet with the surrounding ISM (e.g., 
Brenneman et al.ll2009l ). Our best-fitting continuum model, therefore, has a blackbody and 
only one power-law component at F 
of-fit as the two power-law model. 



1.84^0^0 ci-nd exhibits the same statistical goodness- 



Evidence for excess emission from 2.5 — 7.0 keV was also found by Tiirler et al. , sup- 
porting the presence of a broad Fe Ka line when 3C 273 is in a jet-minimum state. The 
authors quote this excess as significant at the 6a level with an integrated flux of 2.6 ± 
0.4 X 10~'^phcm^^ s^^ , corresponding to an EW = 166 ± 26 eV. These values are consis- 
tent with those repo r ted from previou s observations of the source i n a jet-minimum state 
jKataoka et al.l I2OO2I : IPage et al.l booi IVaqoob fc SerlemitsosI boooh . The breadth of the 
excess is not satisfactorily fitted by a Gaussian line or a diskline component from a non- 
rotating BH. If the line extent is real, the only remaining explanation, according to the 
authors, is that it is emitted from around a near-extreme Kerr BH. The sharp edge of the 
i ron line at 7 keV suggests that the angle of inclination of the accretion disk is 35 — 40° 
(ITiirler et al.ll2006l ). This is contrary to radio observations of the source which have mea- 
sured superluminal velocity in the jet cons istent with it having an inclination angle of < 14° 
to the line of sight (IBallantyne et al.ll2004i ). Such a marked difference between the disk and 
jet inclination angles may imply a warped disk, or perhaps it may mean that the jet incli- 
nation is different on smaller, unresolvable scales. Deeper observations are needed in order 
to resolve this discrepancy. 
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We find evidence for a similar excess around tlie Fe K region in our analysis. Interest- 
ingly, we do not see a narrow emission line at 6.4 keV, but rather the narrow line core seems 
to coincide with the He-like line of Fe Ka at 6.68 ± 0.07 keV with an EW = 49l|^eV. This 
may mean that the gas in this system is significantly ionized, which would not be surprising 
in such an active source. Indeed, this hypothesis is lent credence by fitting a full ref lion 
disk spectrum to the data (which also renders the soft blackbody component obsolete): the 
disk ionization parameter is only constrained to ^ < 5689 erg cm s~^ . Visually, the residuals 
left over after fitting this narrow line suggest the presence of a broad component. Modeling 
this component with our analysis method provides, at best, only a marginal improvement in 
the global goodness-of-fit (Ax^/Az/ = —9/ — 6), perhaps due to the distance of the source 
and the paucity of photons in the data as compared with the other sources in our sample. 
With a 20 ks observation we have collected ~ 3 x 10^ photons from 2.5 — 10 keV. Given more 
observing time and a larger number of (pile-up free) photons, our line analysis would be 
significantly improved. Taking into account these caveats about the robustness of the broad 
line, the kdblur (ref lion) model still provides the best fit, as shown in Tabled though this 
fit is statistically no better than that of a broad Gaussian or laor line. The disk emission 
is highly centrally concentrated with a = 6.O4I206) consistent with the small inner radius 
of emission fit by the model: rmin < 3.73 rg. The X-ray constraints on the inclination angle 
of the disk are consistent with those obtained by Tiirler et al. , within errors: i = 58^22 
degrees. The amount of reflection in this source is difficult to constrain without better pho- 
ton statistics: i^refl = 0.89loi89 • Fig- HI shows the iron line residual, best fit to this residual, 
and relative contributions of the individual model components in the best fit, respectively. 
The diskline and laor fits are shown in Table [TUl while comparisons of the ionized disk 
refiection spectrum with and without relativistic smearing can be found in Table [TTl 



4.2.5. NGC 2992 



NGC 2992 is a Sy 1.9 galaxy {z = 0.0077) that appears to show a broad iron line 
even though it is highly obscured. This source has been the subject of intense study due 
to the variability of its X-ray emission (IGilli et al.l l2000l ). In 1997 and 1998, BeppoSAX 
caught NGC 2992 transitioning from a Compton-thick to a Compton-thin state, resulting 
in an order of magnitude increase in its X-ray luminosity as well as a qualitative difference 
in its spectral appearance, with more disk features (e.g., broad lines) being seen in several 
wavelengths. In the two BeppoSAX pointings, N GC 2992 displays 2 — 10 keV X-ray fluxes 



of 0.63 and 7.4 x 10"" ergcm'^ s"^ , respectively flGiUi et al.ll2000f ) 



We use the May 2003 XMM-Newton observation of NGC 2992, totaling 29 ks, which 
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Fig. 4. — Top left: The 2.5 — 10.0 keV spectrum of 3C 273 fit with a power-law + black- 
body model modified by Galactic photoabsorption. x^l^ = 248/270(0.92). Top right: 
The 2.5 — 10.0 keV best-fit model for 3C 273, including our continuum model and an 
ionized disk reflection spectrum convolved with the kdblur relativistic smearing kernel. 
X^/iy = 229/264(0.86). Bottom: An F(E) plot depicting the relative flux in each of 
the model components in the spectrum. 
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Table 5: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of 3C 273, including 
components and parameter values from 0.6 — 1.5 keV for completeness. As before, the energies 
from 1.5 — 2.5 kcV were not included due to the presence of absorption edges from the 
XMM-Newton mirrors. Error bars are quoted at 90% confidence. We required the zgauss 
line representing the core of the Fe Ka feature to be intrinsically narrow, i.e., a — OkeV. 
Interestingly, this line was found at a moderate ionization state of iron (6.68 keV) rather 
than at the neutral rest-frame energy of 6.4 keV. Redshifts were frozen at the cosmological 
value for the source, in this case 2; = 0.1583. No warm absorption is seen in this object. 
Note the large error bars on the value of R^cPl that we calculate. These are likely an artifact 
of the relatively low photon count for this observation. 
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translates to ~ 3.4 x 10^ photons in the filtered data set. This observation appears to be 
unpubhshed, so all the fit values referenced herein are our own, compiled using the reduction 
and analysis methods detailed in §21 

Upon first inspection, it is immediately clear that NGC 2992 is heavily absorbed below 
~ 2 keV. Though statistical evidence exists for a soft excess, none of the telltale signatures 
of intrinsic warm absorption are present, and the fiux decreases precipitously at soft energies. 
When a simple photoabsorbed power-law model is applied to the continuum, the remaining 
residuals on the soft end indicate that some absorption remains unaccounted for. Including 
a zpcf abs component to model this intrinsic absorption with a partial-coverer in addition 
to the Galactic hydrogen column of iVn = 5.26 x 10^° cm~^ neatly corrects the discrepancy. 
The final continuum parameters are N112 = 5.39^'^o'gg x 10^^ cm~^ with a covering fraction of 
/cov = 95;1^Q% absorbing a power-law of photon index of F = 1.711;q;q4. 

The Fe Ka line is present with a narrow core at 6.43lQQ2keV and EW = 36!l^g^eV. 
Significant residuals remain surrounding the core, however, suggestive of the presence of a 
broad component to the Fe Ka line. Adding a broad Gaussian component at 6.40 keV results 
in a significant improvement in the global goodness-of-fit of Ax^/Az/ = —25/ — 3, and the 
fitted width of o" = 0.20 — 0.51 keV implies an origin of ~ 30 — 190 rg from the black hole 
assuming a Keplerian accretion disk. The fit is only marginally improved by modeling this 
broad component with a diskline, laor or relativistically smeared ionized disk refiection 
model, though the best statistical fit is technically achieved with a kdblur (ref lion) model. 
Nonetheless, due to a lack of significant difference between the rdblur and kdblur fits, 
we cannot distinguish between a spinning and non-spinning black hole in this source, nor 
can we definitively affirm the presence of relativistic effects on the Fe Ka line profile. See 
Table [6] for model parameters and error bars. Fig. [5] shows the iron line residual, best fit to 
this residual, and relative contributions of the individual model components for the best fit. 
The diskline and laor fits are shown in Table [TDl while comparisons of the ionized disk 
refiection spectrum with and without relativistic smearing can be found in Table [TTJ 



4.2.6. NGC 4051 



NGC 4051, like NGC 2992, is a heavily absorbed Seyfert AG N (Sy 1.5; NLSl) seen 



to vary significant ly in fiux over the course of several observations (IGuainazzi et al.l Il996 
Lamer et al.l l2003l ) . The source is at a redshift of z = 0.0023 and has a typical fiux on 



the order of a few times 10~ ergcm" s~ , though as stated above, this fiux can vary 
significantly on a variety of time scales, along with the spectral characteristics of the source. 
Unusually low fiux states in this object can last for weeks to months, during which time the 
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Fig. 5. — Top left: The 2.5 — 10.0 keV spectrum of NGC 2992 fit with a continuum composed 
of a power-law model modified by Galactic photoabsorption as well as intrinsic absorption by 
neutral hydrogen, / v = 1581/1317 (1.20). Top right: The 2.5 — 10.0 keV best-fit model for 
NGC 2992, including our continuum model and an ionized disk reflection spectrum convolved 
with the kdblur relativistic smearing kernel, x^/^ = 1396/1310 (1.07). Bottom: An E^F(E) 
plot depicting the relative flux in each of the model components in the spectrum. 
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Table 6: Best-fitting model parameters for the 2.5—10.0 keV spectrum of NGC 2992, including 
components and parameter values from 0.6 — 1.5 keV for completeness. The energies from 
1.5 — 2.5 keV were not included due to the presence of absorption edges from the XMM- 
Newton mirrors. Error bars are quoted at 90% confidence. Redshifts were frozen at the 
cosmological value for the source, in this case z = 0.0077. No evidence for warm absorption 
is detected. The iron abundance of the ionized disk could not be constrained, so it was held 
fixed at the solar value. 
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X-ray spectrum shows a hard continuum power-law of spectral slope T ^ 1, but is dominated 
by a softer component at lower energies with F ~ 3 (lUttley et al.ll2004j ). A highly broadened 
and redshifted iron line has also been noted in the low flux state of NGC 4051 with RXTE, 
suggesting that reflection features from the accretion disk close to the BH may remain 



const ant in this source in spite of the large variations in continuum properties (jUttley et al. 
2003h . 



NGC 4051 was observed with XMM-Newton in May 2001 for a duration of 117ks, 
yielding ~ 2.5 x 10® photons in the filtered data set. The EPIC-pn results were first reported 
by Mason et al. , and suggested a continuum described by a power-law "pivoting" around 
100 keV, according to a simultaneous observation with RXTE. Ultraviolet emission from 
the Optical Monitor on XMM-Newton was found to lag the X-ray emission by ~ 0.2 days, 
indicating that it is likely reprocessed X-ray emis sion. The X-ray e mission itself showed 
variability on time scales as small as 1 — 2 hours (IMason et al.ll2002l ). These results were 
expanded upon by Pounds et al. , who noted the intermediate flux of the source at this time 
and validated the presence of both a n iron line and a thermal soft excess during the May 
2001 observation (jPounds et al.ll2004l ). Ponti et al. recently re-analyzed this observation as 
well as a lower-fiux pointing from November 2002 and reinforced the veracity of this model, 
while also considering the comp arable efficacy of a model dominated by ionized reflection 
from radii quite close to the BH (jPonti et al.ll2006l ). 



We have also focused on the May 2001 observation of NGC 4051, following the reduction 
and analysis of Mason et al. and Ponti et al. but using modern calibration files. Because we 
do not employ high energy data above lOkeV in our analysis, we find that the underlying 
0.6 — lOkeV continuum is well fit by a single power-law component with photon index F = 
2.36T 

oio- This power-law is modified by Galactic photoabsorption from neutral hydrogen 
with Nil = 1-32 X 10^° cm~^ as well as ionized intrinsic absorption with Nn = 2.77tg5^ ^ 
10^^ cm~^ and ^ = 81^2 erg cm s^^ . A soft excess is seen below ~ 1 keV, as in so many other 
Seyfert galaxies. The best fit to this feature is initially achieved with a thermal blackbody 
component, with kT = 0.10 ± 0.01. 

A narrow 6.40*^;q q3 keV Gaussian was successfully fit to the core of the Fe Ka line in this 
source, though we note significant residuals marking the red and blue wings of a broadened 
line. The narrow core has an equivalent width of EW = 59tii eV. The remaining broad line 
is best fit here with either a diskline or a laor model with EW ~ 35 eV, but in each case 
fitting this single relativistic disk line leaves several parameters poorly constrained. A better 
fit to the entire spectrum is achieved by convolving a static ionized disk refiection spectrum 
model (ref lion) with a relativistic smearing kernel. This more physical model also has the 
advantage of self-consistently accounting for the observed soft excess without the need to 
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include a phenomenological blackbody component. The spectrum of NGC 4051 is equally 
well fit by a kdblur (ref lion) or an rdblur (ref lion) model, though the kdblur model 
is cited as the best fit due to slightly tighter parameter constraints on the inner radius of 
disk emission (rmin = ^■^^-I'^l'i^g)- Full parameter values and 90% confidence error bars are 
listed in Table [71 Fig. [6] shows the iron line residual, best fit to this residual, and relative 
contributions of the individual model components, respectively. The diskline and laor fits 
are shown in Table [TOl while comparisons of the ionized disk reflection spectrum with and 
without relativistic smearing can be found in Table [TTJ 



4.2.7. Ark 120 



Ark 120 is a bright Sy 1 AGN at a redshift of z = 0.0 327, and has a relati vely constant 



2-10 keV X-ray flux of F2-10 ~ 2.50 x 10"^^ erg cm'^ g-i flVauehan et al.ll2004l ). The source 



is radio-quiet, and due to a lack of o bserved evidence f or intrinsic absorptio n. Ark 120 has 
been labeled a "bare" Sy 1 nucleus (jWard et al.lll987l : IVaughan et all 120041 ). Ark 120 has 
been observed by most of the major X-ray obs ervatories. An EXOSAT observation showed 
the source to have a s teep soft X-ray spec trum (ITurner fc Poundslll989l ). as did a subsequent 
i^OiS'^T observation (IBrandt et al.lll993l ). Furthermore, as mentioned above, these X-ray 
observations showed no indicati on of any warm absorption features. Similar find ings were 
seen in ultraviolet observations ( Crenshaw et al.lll999l : ICrenshaw fc Kraemerll200ll ). 



The source was observed by XMM-Newton i n Au gust 2003 for an effective duration of 
looks (~ 2.3 X 10^ photons by IVaughan et al.l (120041 ). We follow the data reduction and 
analysis of these authors, using the latest calibration files, though we restrict our attention to 
the EPIC-pn camera data. Vaughan et al. find that the continuum emission is well described 
by a simple photoabsorbed power-law model with no evidence for complex absorption intrin- 
sic to the system. This finding is substantiated by a comparison to the spectrum of 3C 273 
in its radio-loud state, when the source is dominated by jet emission and displays a very 
fiat, featureless X-ray spectrum lacking intrinsic emission and absorption features. Vaughan 
et al. calculated the ratio between the spectra of Ark 120 and 3C 273, then normalized it 
by a spectral model for 3C 273. The resulting spectrum clearly shows that Ark 120, too, 
possesses a fiat, featureless continuum above ~ IkeV and that its only significant emission 
line is that of Fe Ka at 6.4 keV. The authors use a multiple blackbody component (among 
many other models, with varying degrees of success) for the soft excess to reduce residuals 
left over from the power-law fit below ~ 1 keV, and a slight curvature of the continuum is 
noted at higher energies tha t is thought to indicate the presence of disk refiection in this 
source (I Vaughan et al.l 120041 ). 
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Fig. 6. — Top left: The 2.5 — 10.0 keV spectrum of NGC 4051 fit with a continuum composed 
of a power-law model modified by Galactic photoabsorption, a warm absorber model and 
a blackbody soft excess, x^l^ = 1837/1402(1.31). Top right: The 2.5 - 10.0 keV best- 
fit model for NGC 4051, including our continuum model, a narrow Fe Ka line core and 
an ionized disk reflection spectrum convolved with the kdblur relativistic smearing kernel. 
X^/j^ = 1468/1392 (1.05). Bottom: An F(E) plot depicting the relative flux in each of the 
model components in the spectrum. 
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Table 7: Best-fitting model parameters for the 2.5—10.0 keV spectrum of NGC 4051, including 
components and parameter values from 0.6 — 1.5 keV for completeness. The energies from 
1.5 — 2.5 keV were not included due to the presence of absorption edges from the XMM- 
Newton mirrors. Error bars are quoted at 90% confidence. We required the 6.4 keV zgauss 
line to be intrinsically narrow, i.e., a = OkeV. Redshifts were frozen at the cosniological 
value for the source, in this case z = 0.0023. A warm absorber is present in this object, as 
in many other Sy 1 sources. 
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Our continuum fit is consistent with that of Vaughan et al. , though we find a statistical 
need for only two blackbody components to parameterize the soft excess emission, rather 
than three: kTi = 0.19 ± 0.01 keV and kT2 = O.OSto oo keV. These two temperatures are 
within errors of the first two blackbodies of Vaughan et al. , but we do not detect the third 
component at kT ~ 0.67 keV robustly. The hard X-ray continuum is well represented by 
a power- law of F = 2.111q q7, again consistent with the Vaughan et al. value within errors. 
The value of the Galactic absorbing column is taken as A^h = 1-26 x 10^^ cm~^ . 

Vaughan et al. detected the presence of broad and narrow components to the 6.4 keV 
Fe Ka line (-EPVbroad ~ 100 eV, -EPVnarrow ~ 40 eV). The narrow component could easily be 
fit with a 6.4 keV Gaussian profile, as was our Fe Ka narrow core, but Vaughan et al. found 
that a diskline model worked best for their broad component (rmin ~ 144 Tg). Interestingly, 
this broad component fit to a rest-frame energy of 6.56 keV ra ther than 6.4 keV, w ith the 



latter yielding a substantially worse global fit (Ax^ = +21.5) (IVaughan et al.l 120041 ). This 
energy corresponds to the line emission being dominated by intermediately ionized iron 
(Fexx-xxil), me aning that resonant trapping and the Auger effect should destroy the line 



flRoss et al.lll996h 



In an effort to solve this puzzle we have fit our iron line profile with three narrow 
Gaussian components, the first two at 6.40 ±0.01 keV and T.OOlg'os keV, representing neutral 
iron and a likely blend of highly ionized H-like Fe Ka and Fe K/3, respectively. The lines have 
-E'W^coid = 36l8^eV and EWio^ = 28lgeV. After an initial attempt to fit the broad residual 
remaining, there is evidence for a third narrow line as well at 6.751^qq7 keV {EW = lA^^eV), 
indicating the presence of moderately ionized iron. The addition of this component improves 
the fit by Ax^ = —19 for two additional degrees of freedom. This fit leaves a residual 
broad feature nicely centered at 6.4 keV, which renders this scenario physically consistent 
with other systems and removes the need to explain the presence of the broad line of iron 
in an intermediate state of ionization. We then achieve a best fit to the broad component 
using a Schwarzschild rdblur (ref lion) model, which also nicely accounts for the additional 
curvature noted in the spectrum at higher energies by Vaughan et al. , as well as the soft 
excess we had initially modeled with phenomenological blackbody components. The best-fit 
parameter values and error bars are presented in Tabled Fig. [7] shows the iron line residual, 
best fit to this residual, and relative contributions of the individual model components for the 
best fit. The diskline and laor fits are shown in Table [T0| while comparisons of the ionized 
disk reflection spectrum with and without relativistic smearing can be found in Table [TTl 
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Fig. 7. — Top left: The 2.5 — 10.0 keV spectrum of Ark 120 fit with a continuum composed 
of a power-law model modified by Galactic photoabsorption and a two-zone blackbody soft 
excess, x^l^ = 2153/1428 (1.51). Top right: The 2.5 - 10.0 keV best-fit model for Ark 120, 
including our continuum model and an ionized disk reflection spectrum convolved with the 
rdblur relativistic smearing kernel, x^/v = 1511/1418(1.07). Bottom: An E^F(E) plot 
depicting the relative flux in each of the model components in the spectrum. 
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Table 8: Best-fltting model parameters for the 2.5 — 10.0 keV spectrum of Ark 120, including 
components and parameter values from 0.6 — 1.5 keV for completeness. The energies from 
1.5 — 2.5 keV were not included due to the presence of absorption edges from the XMM- 
Newton mirrors. Error bars are quoted at 90% confidence. We required the zgauss lines to 
be intrinsically narrow, i.e., a = OkeV. Redshifts were frozen at the cosmological value for 
the source, in this case z = 0.0327. 
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4.2.8. Fairall 9 



Fairall 9 is a radio-quiet Sy 1 galaxy with an elliptical companion, both at a moderate 
redshift of z = 0.047. This source has not been observed to undergo very large changes 
in X-ray flux over the time it has been o bserved: the typical 2 — lOkeV flux of Fairall 9 
is F2-10 ~ 1.5 — 5.0 X 10^^^ erg cm"^ s~^ (IReynoldsl 119971 : iGondoin et al.ll200ll ). Previous 
X-ray observations with ASCA show that the source has a continuum well- described by a 
photoabsorbed power-law with a fluorescent Fe Ka line and a high-energy tail. The latter 
is a common signature of disk reflection, as has been discussed previously w ith respect to a 



number of other sources in this Section (IRevnolds 



1997 



Nandra et al.l 119971 ). A soft excess 



has also been detected below ~ 2keV JPounds et al.lll99iL_ though no strong evidence of 



warm absorption has been observed in the soft spectrum (IGondoin et al.ll200ll ). 



Th e source was observ ed with XMM-Newton in July 2000 for an effective duration of 
29 ks (jjansen et al.ll200l[ ). resulting in a flltered EPIC-pn data set with ~ 2. 0x10^ photons. 



The E PIC-pn data were presented in full by Gondoin et al. the following year (IGondoin et al. 



2OOII ). We have followed the general procedures for data reduction and analysis discussed 
by Gondoin et al. , using the most up-to-date calibration flies and software as we have for 
all other sources analyzed in this work. 



Gondoin et al. flnd that continuum of Fairall 9 is best modeled with a photoabsorbed 
power-law (P = 1.791;q ;[2)5 which the absorbing column density of neutral hydrogen gas 
is approximately equal to its Galactic value along the line of sight to this source (iVn = 
3.19 X 10^° cm~^ ). A soft excess is noted below ~ 2 keV, and an RGS analysis yielded a best- 
fitting blackbody model for this component: kT = 0.17!tQQ7keV. Incorporating reflection 
from an ionized disk (pexriv) with solar abundances and a reflection fraction of i?refl = 1-0, 
the authors flnd that the high-energ y tail is quite well modeled by a disk at an inclination 
angle of z = 26° (IGondoin et al.ll200ll ). We have performed a slightly different flt to the EPIC- 
pn data only, replacing the pexriv component with a relativistically-blurred ref lion model, 
which also accounts naturally for the soft excess without employing a phenomenological 
blackbody component. The column density of neutral hydrogen absorption was again frozen 



Itoil, and the 



at the Galactic value of = 3.19 x 10^" cm"^. We found that P = 1. 
disk inclination angle flt to z ~ 54°, though its value could not be properly constrained. 
Unfortunately, due to the relatively small number of photons in the 0.5 — 10.0 keV spectrum, 
parameter constraints were more problematic for this source than for many of the others in 
our sample. 

Gondoin et al. also note the presence of a ~ 6.4 keV Fe Ka line in the spectrum with 
EW ~ 120 eV, which they modeled with a Gaussian component. Due to its relatively narrow 
proflle, the line is suggested to originate from low-ionization material orbiting relatively far 
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out in the disk. Additionally, an Fe K absorption edge is also noted at 7.64 keV with 
r = 0.18, consistent with reflection from cold, optically thick material. We find that the 
inclusion of such an edge is not statistically robust in our model once reflection is included 
via ref lion (which incorporates the iron edge), but we do find the narrow 6.40lQ'Q2keV line 
in emission and model it with a narrow Gaussian feature representing the core of the Fe 
Ka line. This feature has an EW = 115l2oeV, but even after its inclusion residual features 
remain indicating the presence of other unmodeled narrow lines as well as a potential broader 
iron line component as evidenced by a broadened red wing. Adding in a second narrow line 
at E = 7Mtom 

keV {EW = 20^QeV), representing a contribution from Fe K/3 emission 
(and so a flux of 0.16x that of Ka), resulted in a marginal statistical improvement in fit 
(Ax^/Az/ = —8/— 2), though the visual improvement in fit was readily apparent. Attempting 
to model the remaining residual feature with a single broad line produced unsatisfactory 
results, however. As such, we added a third narrow Gaussian line at 6.78lQ;Q5keV {EW = 
29^*^9^ eV), representing ionized iron, which further improved the fit by A^^/Ai/ = —7/ — 2. 
While both the ionized iron line and Fe K/3 line did not result in statistically significant 
improvements in fit, according to the F-test, they nonetheless produced substantial visual 
improvements to the fit and nicely modeled the residual features that could not be accounted 
for by a broad Fe Ka line component. A residual emission line-like feature remains around 
~ 7.1 keV in the observed frame, which likely represents a line from Nix-xvi, but adding in 
an additional Gaussian to model it results in no net improvement in fit so it is left unmodeled. 

It should be noted that Gondoin et al. did not report the detection of a broad Fe Ka line 
in their analysis, and that our attempts to include one show that, while the inclusion of a 
broad component at 6.40 keV does improve our fit by A^^/Az/ = —13/ — 2, a relativistic disk 
line model such as diskline or laor shows no improvement over a simple Gaussian with 
a = 0.16 ± 0.11 keV. Further, the parameter values of the diskline and laor components 
are largely unconstrained, as one might expect. Even so, our best statistical fit was achieved 
with a rdblur (ref lion) model, albeit with the emissivity index and inclination angle of 
the disk unconstrained. The inner radius of disk emission was constrained to rmin < 29 r^. 

Best-fit values and error bars for all the parameters in this model are listed in Table [91 
Fig. [H] shows the iron line residual, best fit to this residual, and relative contributions of 
the individual model components, respectively. The diskline and laor fits are shown 
in Table [TOl while comparisons of the ionized disk reflection spectrum with and without 
relativistic smearing can be found in Table [TTl 
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Fig. 8. — Top left: The 2.5 — 10.0 keV spectrum of Fairall 9 fit with a continuum composed 
of a power-law model modified by Galactic photoabsorption and a blackbody soft excess. 
X^/z/ = 816/699 (1.17). Top right: The 2.5 — 10.0 keV best-fit model for Fairall 9, including 
our continuum model and an ionized disk reflection spectrum convolved with the rdblur 
relativistic smearing kernel. xV^ = 679/691 (0.98). Bottom: An Y?Y{Yj) plot depicting the 
relative flux in each of the model components in the spectrum. 
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Table 9: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of Fairall 9, including 
components and parameter values from 0.6 — 1.5 kcV for completeness. The energies from 
1.5 — 2.5 keV were not included due to the presence of absorption edges from the XMM- 
Newton mirrors. Error bars are quoted at 90% confidence. We required the 6.4 keV zgauss 
line to be intrinsically narrow, i.e., a — OkeV. Redshifts were frozen at the cosmological 
value for the source, in this case z — 0.0470. A warm absorber is present in this object, as 
in many other Sy 1 sources. 
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5. Results 

We have compiled the relevant parameter constraints from all of our spectral fits to eight 
Sy 1 sources: MCG-5-23-16, Mrk 766, NGC 3783, NGC 4051, Ark 120, Fairall 9, NGC 2992, 
and 3C 273. Table [10] presents the results of fitting the broad iron line with a simple 
diskline or laor model as a first approximation. Table [TT] displays the results of fitting 
a more detailed, physical model consisting of many reflection features from the accretion 
disk as a whole, both with and without relativistic smearing via rdblur and kdblur. Note 
that, for nearly every source, a model including reflection from an ionized disk convolved 
with a relativistic smearing kernel provides the best statistical fit, as judged by the change 
in chi-squared per additional degree of freedom in the model (i.e., the F-test), tightness of 
the parameter constraints and reduction in the residuals from the data/model ratio. 

Previous studies have indicated that broad iron lines may be present in up to 42% of 
AGN observed with the XMM-Newton/FjPlC-pn camera that have > lO'* counts (IGuainazzi et al 



20061 ). To assess whether evidence exists for a relativistically broadened iron line in each 
of our data sets, we first considered the diskline and laor fits in which BH spin is fixed 
at a = 0.0 and a = 0.998, respectively. While each source demonstrated a significant im- 
provement in its global fit as compared to a model fitting only a narrow iron line core with 
a Gaussian component, the most important parameter to evaluate in this case is the inner 
radius of disk emission. Roughly speaking, as stated in §3.2^ if r^^^ < 20 rg we can say 
with some confidence that there is substantial emission from the inner part of the accretion 
disk where relativistic effects (such as BH spin) become important in shaping the overall 
iron line profile. Out of our eight sources, Fairall 9 and Ark 120 had very poor constraints 
on rmin based on their diskline and laor fits. Of the six sources where r^^^ could be 
constrained, only MCG-5-23-16 did not show a substantial contribution from r^in < 20 rg, 
either in the diskline or laor fit. All of our sources showed little difference between the 
diskline and laor model fitting results, statistically, as judged by their reduced values. 
The four sources that did show very modest reductions in reduced with a laor model 
were MCG-5-23-16, Mrk 766, NGC 2992 and 3C 273. Of these, MCG-5-23-16 showed the 
most significant improvement in fit with laor as opposed to diskline, though it is doubtful 
that this difference represents evidence for a spinning BH, given the constraint on the inner 
disk radius in this source. 

We have also examined the question of how robust the presence of a broad line is in the 
spectrum when reflection is included. Beginning with a ba se continuum model including a 



static ionized disk reflection spectrum (IRoss fc Fabianll2005l ). we noted the overall goodness 



of-fit and the residual features remaining, especially around the iron line region. We then 
convolved this model with relativistic effects, first using the rdblur smearing kernel, then 
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substituting in kdblur (Table [3]). Due to the nature of the BH spin parameter space, the 
difference between an a = 0.8 and an a = 0.998 BH is much more pronounced than the 
difference between an a = 0.0 and an a = 0.8 BH. Therefore, while a significant improvement 
in the global fit with either model indicates that relativistic smearing is important in the 
system, further significant improvement with kdblur instead of rdblur suggests that not 
only is relativity important, but in the best-fit scenario the BH spin likely tends toward the 
maximum theoretical value. 

The rdblur and kdblur fits were at least as good or better than the diskline and 
laor model fits in each source. The difference is greatest in the cases of MCG-5-23-16 
and NGC 4051. These smearing models also returned similar results in terms of parameter 
constraints and overall fit to the iron line profile in each source, and in nearly all cases also 
produced an improvement in the overall goodness-of-fit over the non-smeared case. More 
marginal improvements were seen in Fairall 9, NGC 2992 and 3C 273; not coincident ally, 
these are the three sources for which we have the lowest photon count (each has < 4 x 10^ 
photons from 2 — lOkeV). Note that it was also difficult to constrain several other parameters 
in these AGN: the emissivity index and inner radius of disk emission, particularly, and in 
the case of Fairall 9 even the inclination angle of the disk to our line of sight. Normally 
inclination is well-constrained by the position of the blue wing of the Fe Ka line profile, but 
Fairall 9 has the fewest counts of any of our sources, and so the blue wing cannot easily 
be separated from the underlying continuum. This underscores the importance of photon 
statistics in evaluating the relevance of relativistic effects in an AGN spectrum. Out of the 
five sources that showed significant improvement with relativistic smearing, none showed 
a significant statistical preference for either rdblur or kdblur according to their reduced 
chi-squared values, implying that we are not able to robustly distinguish between a non- 
spinning and a rapidly-spinning BH in any of our sample members. MCG-5-23-16, Ark 120, 
NGC 2992 and 3C 273 all show a modest preference for kdblur over rdblur, while NGC 3783 
shows the opposite. The other sources show equally good fits to each model. 



We note that the reflion model (iRoss &: Fabian! |2005| ) does not include a formal pa- 



rameter for the reflection fraction of the disk, so in order to facilitate comparison between our 
results and those from s tudie s using other models to parameterize reflection (e.g., pexrav. 



Magdziarz fc Zdziarskil (Il995l )). we have estimated R^ea in all of our models shown in Ta- 
ble [11] and in §4. To make these estimates, we have assumed that for i^rcfi = 1, half of the 
radiation from the hard X-ray power-law is received and reflected by the disk, with the other 
half traveling directly to the observer. Reflection fractions of greater than unity then indicate 
that more than half of the power-law radiation illuminates the disk, implying anisotropic 
emission from the power-law source or perhaps strong light-bending of this emission near 
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the BH. Operationally, we compute RreH using the following formula: 



iVpo Fpo V30, 

Here Nj-cQ and iVpo denote the normalizations of the ref lion and power-law components 
from our best-fit spectral model, (^) is the best-fit ionization parameter characterizing the 
disk refiection, and F^cH and Fpo are the total fiuxes contained in the ref lion and power-law 
components for unity normalization over the full wavelength range (0.001 — 1000 keV). These 
reflection fractions are often poorly constrained due to the lack of simultaneous high energy 
data above lOkeV, which would enable us to constrain the amount of reflection (as well as 
the iron abundance and ionization of the disk) more effectively by measuring the Compton 
hump, which typically peaks at 20 — 30keV. Having this high energy data would therefore 
allow us better parameterize many properties of the disk, enabling better constraints on 
the nature of the spacetime surrounding the BH (e.g., rmin and the BH spin). For this 
reason, part 11 of this work will utilize Suzaku/ XIS+B.XD data to establish robust BH spin 
constraints in these AGN. 
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Table 10: Comparison of the spectral fitting results for the 2.5 — 10.0 keV spectra using the diskline and laor models. 
For each source the top row represents the diskline range in parameter value while the bottom row represents the 
laor range in parameter value. Ranges are given at 90% confidence. A dashed line rather than a numerical value 
indicates that the parameter could not be statistically constrained in a formal error analysis. 
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Table 11: Comparison of the spectral fitting results for the 2.5 — 10.0 keV spectra using the ionized reflection model 
reflion, rdblur (ref lion) and kdblur(ref lion). The latter two models convolve the reflection spectrum with 
smearing kernels from a non-spinning and maximally-spinning BH, respectively. For each column of data, the top row 
represents the reflion range in parameter values, the middle row shows the rdblur (ref lion) values and the bottom 
row represents the kdblur (ref lion) values. For the reflion component, the value of the incident power-law spectral 
index is set equal to that of the hard x-ray source. Reflection fraction is estimated as described above in the text, based 
on the relative fluxes and normalizations of the power-law and reflion components of this best-fit model. Parameter 
value ranges are given at 90% confidence. 
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6. Discussion and Conclusions 

This paper explores the use of broad iron emission lines to constrain the importance 
of relativistic effects in a sample of AGN. Following our previous work on MCG-6-30-15 
(BR06), this study of an additional eight prominent Sy 1 AGN observed to have prominent 
iron lines has provided us with some intriguing data on this subject and has also put us in 
an excellent position to apply our spectral fitting techniques to many other sources in the 
future. 

Though this sample is far too small to allow us to draw any robust conclusions about 
the nature of the spacetime in the nuclei of Sy 1 galaxies as a whole, this work has allowed 
us to design a viable method for evaluating the presence of a broad iron line in the data 
and using a sequence of models to assess the importance of relativistic smearing in shaping 
the spectral features emitted from the accretion disk around the BH. We have determined 
that four out of our eight sources show improvement in their goodness-of-fit parameters 
when relativistic effects from the inner disk are taken into account (one more shows at least 
marginal improvement). Of these four, two show more improvement with a smearing model 
designed for a maximally-rotating BH than for a non-spinning BH, which may suggest that 
BH spin is non-zero in these objects. 

Access to energies greater than the XMM limit of lOkeV will help to better constrain 
the disk parameters these sources, in turn allowing for constraints to ultimately be placed 
on BH spin. This endeavor will be the focus of our follow-up Suzaku work on these AGN. 

Based on our XMM data analysis we have arrived at three important conclusions that 
we hope will inform subsequent broad iron line and BH spin surveys: 



Guainazzi et al. ( 20061 ) found that the fraction of AGN displaying broad iron lines in 



XMM-Newton/'EPIC data increases dramatically when the spectra have more than 
~ 10^ counts. We find that EPIC spectra with ~ 10^ counts are required in order to 
use these iron lines to assess the importance of relativistic effects in these sources. 

Five out of nine sources from our broad Fe K line sample (including previous work on 
MCG-6-30-15) showed a significant improvement in their global goodness-of-fit when 
relativistic smearing was added to a static ionized disk reflection spectrum. 

Though the XMM-A'^ewton/EPIC-pn instrument allows us to resolve the broad iron 
line in excellent detail, its lOkeV energy cutoff limits our ability to accurately model 
the disk reflection continuum due to its inability to measure the peak energy of the 
Compton reflection hump at ~ 20 — 30 keV. 
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Purther progress in this direction needs more data with newer, more sensitive, broader 
bandwidth instruments in order to provide constraints on BH spin. Suzaku, with its enhanced 
spectral resolution in the XIS detector and the HXD/PIN detector simultaneously covering 
the energy range from 10 — 60keV, provides great advantages for this type of study. Part 
II of this work (forthcoming) will analyze archival Suzaku observations of these sources, 
enabling us to measure the peak of the Compton hump and break the degeneracies between 
the parameters of the disk and those governing the spacetime near the BH. Suzaku data will 
therefore allow us to constrain BH spin using our kerrdisk and kerrconv models. 
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